Magnetostriction-driven multiferroicity of MnTe and MnTe/ZnTe epitaxial films 
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Here we demonstrate that MnTe epitaxial films with zinc-blend structure and MnTe/ZnTe multi- 
layers should show ferroelectric polarization in antiferromagnetically (AFM) ordered state and thus 
belong to multiferroics. Spontaneous ferroelectric polarization results from the bending of highly 
ionic Mn-Te-Mn bonds induced by magnetostrictive shear strain. Orientation of ferroelectric polar- 
ization is coupled with orientation of AFM vector and thus can be controlled by application of the 
external magnetic field. Due to the clamping of electric and magnetic order parameters, domain 
structure in MnTe is governed by two mechanisms: depolarizing field produced by electric dipoles 
and destressing field produced by magnetoelastic dipoles. The values of monodomainization electric 
and magnetic fields depend upon the sample shape and diminish with the film thickness. Magne- 
toelectric nature of the domains make it possible to visualize the domain structure by linear and 
nonlinear optical methods (Kerr effect, second harmonic generation technique). 

PACS numbers: 75.50.Ee Antiferromagnetics 75.50.Pp Magnetic semiconductors 75.80.-|-q Magnetomechan- 
ical and magnetoelectric effects, magnetostriction 77.80.-e Ferroelectricity and antiferroelectricity 
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I. INTRODUCTION 

Multiferroics are a class of a single phase or compos- 
ite materials with coexisting magnetic and ferroelectric 
ordering. High susceptibility of multiferroics to both elec- 
tric and magnetic fields which enables electrical control 
of magnetic state and vice versa makes them interesting 
from fundamental and applied points of view. 

Some of multiferroics are ferroelectric and ferromag- 
netic (see, e.g., Refs. [111111) and could be used as a storage 
media with high information capacity or as functional- 
ized materials for electronic devices. Others (mainly rare 
earth compounds with Mn ions like RMnOa or RMn205, 
where R — Tb, Dy, Ho, Y, etc.) show ferroelectric and 
antiferromagnetic ordering^i^iii^ and may find an appli- 
cation as mediators for indirect electric control of ferro- 
magnetic state of an adjacent layer.— 

Cross-coupling between the magnetic and ferroelectric 
properties may be caused by different mechanisms. In 
type-I multiferroics (according to classification proposed 
in the recent review^*^) ferroelectricity and magnetism 
have different origin and coexist only in a certain tem- 
perature range. Materials classified as type-H show fer- 
roelectricity only in a magnetically ordered state. Until 
now such a combined ferroelectric-magnetic ordering was 
observed in crystals which posses at least small noncom- 
pensated magnetization, i.e. in ferro-. ferri- and weak 
ferromagnets (see, for example, Refs. l^ fTllll^llsl . etc.). 

Two microscopic mechanisms of the magnetically- 
induced ferroelectricity originate from the exchange 
interactions^, namely, i) anisotropic exchang o^^d^ 
(Dzyaloshinskii-Moria interaction, DMI) and ii) exchange 
striction^^ (magnetoelastic coupling to the lattice). In 
the first case antisymmetric DMI activated by non- 
coUinear spin ordering breaks inversion symmetry and 
induces electric polarization through the concomitant lat- 
tice and electronic distortion. Some typical examples of 
the compounds which show DMI are listed in Table [H 



The second mechanism may play the dominant role 
in the antiferromagnets with superexchange interactions 
between the magnetic ions mediated by bridging via 
nonmagnetic anions, like DyMn205 and YMn205, see 
Table IT] In this case the value of exchange integral 
strongly depends upon an angle between anion-cation 
bonds. In the crystals with the competing antiferromag- 
netic (AFM) interactions even small bond-bending may 
reduce the exchange energy, stabilizes long-range AFM 
ordering and produce nonzero electric polarization. 

In the present paper we predict multiferroicity caused 
by magnetoelastic mechanism in the epitaxially grown 
MnTe/ZnTe films and heterostructures that belong to 
the family of II- VI semiconductors^^ well suited to many 
optoelectronic applications in the infrared and visible 
rangei^. Our hypothesis is based on the following facts. 

1. Epitaxial layers of MnTe, ZnTe, CdTe and corre- 
sponding heterostructures posses zinc-blend (ZB) 
structure consistent with piezoelectric activity. In 
particular, isomorphous to MnTe nonmagnetic 
films of ZnTe, CdTe, ZnTe/CdTe(lll) can produce 
a macroscopic electric polarization when stressed 
or strained^S. 

2. Predominant mechanism of the magnetic interac- 
tions between Mn^+ ions is superexchange via Te 
anions^!^i2^ that favours antiferromagnetic order- 
ing and noticeably varies with Mn-Te-Mn bond 
bending. 

3. Magnetoelastic coupling in MnTe is rather strong 
as can be deduced from step-wise variation of lat- 
tice parameters at the Neel temperaturej^i^i 

4. MnTe is a wide-gap (band gap is 3.2 eV, as re- 
ported in Ref. (iSj) semiconductor with vanish- 
ingly small concentration of the mobile carriers at 
low temperaturesi^ Thus, ferroelectric polariza- 
tion should not be seriously affected by screening by 
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the mobile charge. In MnTe/ZnTe heterostructures 
such a screening can be avoided by a proper choice 
of the superlattice parameters and/or by doping. 

We argue that the value of spontaneous ferroelectric 
polarization in MnTe can be as large as 60 nC/cm^ 
and thus is comparable with polarization of many other 
antiferromagnets (see Table HI. Moreover, due to the 
direct coupling between magnetic and ferroelectric or- 
dering, one should expect strong magnetoelectric effects 
(i.e., induction of polarization by a magnetic field) in this 
material. Magnetoelastic mechanism is responsible also 
for the nonlinear optical effects and opens a possibility 
to visualize antiferromagnetic (and corresponding ferro- 
electric) domains by use of second harmonic generation 
(SHG) technique (see, e.g, Refs. [UliS). 

The paper is organized as follows. In SecHIlwe describe 
crystal and magnetic structure of MnTe/ZnTe films and 
in Sec lIIII give some intuitive reasons for appearance of 
spontaneous electric polarization below the Neel point. 
In Sec lIVI we calculate the value of spontaneous polariza- 
tion on the basis of phenomenological model. Sec|V]is de- 
voted to the discussion of AFM and ferroelectric domains 
and methods of the domain structure control. In Sec lVIl 
we discuss possible magnetooptical effects that could be 
helpful in visualization of the magnetic and ferroelectric 
structure of MnTe. In the last Sec. the summary of the 
results obtained is given. 

II. STRUCTURE AND MAGNETIC 
PROPERTIES OF MNTE 

Bulk-like films of paramagnetic MnTe grown by MBE 
technique have zinc-blend structure?- shown in Fig. [T] 
(symmetry group is F43m) . The d — d exchange interac- 
tion between the nearest localized Mn^+ spins {S = 5/2) 
is accomplished via Mn-Te-Mn bonds^i and thus IS an- 
tiferromagnetic. In the strain-free fee lattice AFM ex- 
change between the nearest neighbors (NN) is frustrated. 
Nevertheless, neutron diffraction experiments^^ii^ reveal 
appearance of the long-range AFM type-Ill structure be- 
low the Neel temperature (Tjv — 65 K). As shown in 
Fig. [2l this spin arrangement consists of AFM sheets par- 
allel to (001) plane (represented in the plot by the shaded 
parallelograms). Stabilization of long-range ordering is 
attributed to the influence of weak NNN coupling2^ and 
to misfit-induced strain of the magnetic layer— which 
in the case of MnTe/ZnTe corresponds to elongation^^ 
(c/a = 1.06) in the direction of film growth. 

It should be noted that in the described AFM-III 
structure the coupling energy between the spins in the 
adjacent (001) sheets sums to zero, so, mutual rota- 
tion of spins in different layers does not change the 
energy of exchange interactions^. Stabilization of the 
collinear AFM-III structure can be due to the pro- 
nounced mag 

netostriction^?, (0.3%) that removes degen- 
eracy between [100] and [010] directions within the film 
plane. 




Figure 1: (Color online) Structure of zinc-blend MnTe. Mn 
- magenta (large) spheres, Te - blue (small) spheres. Inset 
(small cube) shows orientation of Mn-Te bonds. 



[001] 




Figure 2; (Color online) Collinear AFM-III structure in the 
tetragonally distorted fee spin lattice. Film growth direction 
is parallel to [001] crystallographic axis. Mn^+ ions with op- 
posite spin direction are shown in different colours (opposite 
arrows) . 



It should be also mentioned that the symmetry of AFM 
state (corresponding group is generated by the rotation 
2 [001] and translation [i, i,0] both combined with time 
inversion 1') allows the existence of macroscopic electric 
polarization vector oriented along the [001] axis and for- 
bids existence of macrosopic magnetization. 
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Table I: Values of spontaneous electric polarization, p(=P°")j Neel temperature, Tjv, type of AFM ordering (MO), microscopic 
mechanism and type (according to classification Ref Hol) of ferroelectricity for some multiferroics (MF). 



Compound 


p(spon) (nC/cm^) 


Tat (K) 


Type of MO 


Mechanism 


Type of MF 


Source 


BiFeOs 


12-10" 


673 


cycloid 


octahedra rot. 


I 


Ref.9 


BiSrFeOg 


lO-lO" 


643 


noncoUinear 


DMI 


I 


Ref.28 


DyMn205 


150 


39 


noncoUinear 


ME 


II 


Ref.29^ 


YMnaOs 


100 


45 


coUinear 


ME 


II 


Ref.rz 


MnTe (ZB) 


60 


65 


coUinear 


ME 


II 


this work 


TbMnOs 


50 


27 


spiral 


DMI 


II 


Ref.31 


TbMnaOs 


40 


37 


noncoUinear 


DMI 


II 


Ref.^ 


NisVaOs 


12.5 


6.3 


noncoUinear 


DMI 


II 


Ref.32^ 


MnW04 


4 


12.7 


noncoUinear 


DMI 


II 


Ref.I 


CasCoMnOe 




13 


coUinear, Ising 


ME 


II 


Ref.34 



III. MULTIFERROICITY INDUCED BY 
MAGNETOSTRICTION: INTUITIVE 
CONSIDERATIONS 

The role of magnetoelastic coupling in the formation of 
AFM-III structure can be traced from the following qual- 
itative considerations. It was already mentioned that the 
value of Mn-Te-Mn angle (f> is the key quantity in deter- 
mining the exchange coupling constant Jd-d between the 
NN Mn^+ ions: upon decreasing the (j) the strength of 
AFM interaction increases. For small deflection from an 
ideal 00 = 109.5° angle peculiar to fee lattice, this depen- 
dence can be approximated as follows (from the results 
of Bruno and Lascaray'^I) 



Jd-dW^~3A5 + 0.lS5{(j)-M, K. 



(1) 



Now, let us turn to Fig. [3] which illustrates the effect of 
tetragonal strain induced by the presence of ZnTe layers. 
In a nondeformed cubic lattice all the bonds make an 
"ideal" angle (f)Q and all NN interactions are equivalent. 
Tetragonal distortion in [001] direction removes degener- 
acy between in-plane (001) and interplane NN exchange. 
In the case of elongation the angles for pairs 12, 34 are 
smaller than those for the other pairs: 



f34 



< 



'13 



(2) 



Typical values for in-plane, (/)in, and interplane, (t>ont, an- 
gles calculated from geometrical considerations for some 
MnTe/ZnTe super lattices are listed in Table [Til In the 
last two columns of the Table we give the values of the ex- 
change constants calculated from Eq.(IT]). It can be easily 
seen that NN exchange interaction favors AFM ordering 
for the atoms within (001) plane. At the same time, the 
bonds between the atoms in neighboring planes are still 
equivalent and so are frustrated. 

Further deformation within (001) plane (see Fig. [4] a, 
b) leads to optimization of the exchange in the interplane 
FM and AFM bonds by strengthening exchange interac- 
tions with the "right" sign and weakening those with the 





(b) 



Figure 3: (Color online) Effect of the mismatch-induced 
stresses. Blue (1, 3) and magenta (2, 4) spheres correspond 
to Mn atoms with opposite directions of magnetization, (a) 
Nondeformed (cubic) cell, atoms in 1-4 positions are equiv- 
alent and angles between all the Mn-Te bonds are the same 
and equal 109.5°. (b) Elongation in [001] direction removes 
degeneracy between angles. 



"wrong" sign. Namely, elongation in [100] direction (cor- 
responding strain component u^x — Uyy > 0) results in 
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(a) 



9ii 



(b) 




Figure 4: (Color online) Magnetostriction in (001) plane. 
Pairs of antiferromagnetically coupled Mn atoms 1, 2 (hollow) 
and 3,4 (filled) belong to different atomic planes. (a) Nonde- 
formed state, angles between Mn-Te-Mn bonds take only two 
different values: (jiin (pairs 12 and 34) and (^out (pairs 13, 14, 
23 and 24). (b) Elongation in [100] direction removes degen- 
eracy between the angles 0i3 — 024 > 0out > 023 ~ 4>i4,- (c) 
Shear strain Uxy removes degeneracy between the angles 0i2 
and 034. 



the following difference between the exchange constants 

0.27axy 



Jd~d{(kii) - Jd-d{4>i'i) = 



^yy 



), K, 



(3) 

where a^y and c are in-plane and interplane lattice pa- 
rameters, correspondingly. 

Thus, magnetostriction can stabilize an AFM-III struc- 
ture even in approximation of NN exchange interactions. 
On the other hand, magnetostriction plays an important 
role in formation of ferroelectric ordering, as we will show 
below. 

In the case of strong spin-orbit coupling some compo- 
nents of magnetostrictive tensor could depend upon the 
mutual orientation of a localized magnetic moment and 
a direction of anion-cation bond. This situation in appli- 
cation to MnTe is illustrated in Figs. Hj;, [5] Shear strain 
Uxy in (001) plane makes AFM bonds between 1-2 and 
3-4 pair inequivalent: 



\/2axy 3c^ 



^34 



2a% + c2 



(4) 



The displacement A£ between the positive (Mn^+) and 
negative (Te^~) ions is thereby equal to 



9^ 



''xy ■ 



(5) 



For the totally ionic bond (with the effective charge 2e) 
and spontaneous strain u^y = 3 • 10~^ such a displace- 
ment induces a local dipole moment cx 0.12 D. Accu- 
rate calculations of ferroelectric polarization induced by 
spontaneous magnetostriction will be given in the next 
SecHVl 





Figure 5: (Color online) Two types of AFM (and ferroelec- 
tric) domains, (a), (c) Mutual orientation of the sublattice 
magnetization (vectors Mj, j = 1,4) and shift vector (uj, 
j = 1,4) of Mn atom, (b), (d) Arrangement of atoms after 
the shear strain. Nonequivalence of 1-2 and 3-4 bonds results 
in appearance of ferroelectric polarization P. 

Figinireveals one interesting feature of the tetragonally 
distorted MnTe system: AFM structure can be imple- 
mented in two types of domains that have different (per- 
pendicular) orientation of local magnetic moments (with 
respect to crystal axes) and opposite direction of electric 
polarization. AFM domains are clamped with ferroelec- 
tric ones. This means that macroscopic state of a sample 
can be controlled by application of either magnetic or 
electric field (or both) . Behavior of the domain structure 
in the presence of the external fields will be discussed in 
Sec. El 



IV. POLARIZATION INDUCED BY AFM 
ORDERING: MODEL AND CALCULATIONS 

Formally, AFM type-Ill structure can be described by 
the following distribution of magnetization vectors M(rj) 
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Table 11: Angles between Mn-Te-Mn bonds for in-plane NN (i^in) and interplane NN (0out) and corresponding values of exchange 
integrals Jd-d(0in), Jd-d{4'ovLt) (K) calculated from geometrical considerations for different superlattices. In-plane, a^y, and 
interplane, c, lattice parameters (A) are taken from the experimental paper Ref|23|. 



Type Oxy C 0out Jd-d{<pin) Jd-d{(f>out) 

MnTelO/ZnTelS 6.130 6.570 105.7 111.4 -3.20 -3.97 

MnTe20/ZnTel8 6.183 6.505 106.7 110.9 -3.27 -3.84 

MnTel30/ZnTe330 6.150 6.470 106.7 110.9 -3.28 -3.83 

MnTe 6.346 6.346 109.5 109.5 -3.45 -3.45 



at Mn sites: 

M(r,) =L(q)exp(zqr^-)- (6) 

Three equivalent orientations of the structure vector 
q = (27r/a, 0, 7r/a) (where a is a lattice constant), and 
two orientations of L vector (parallel to [100] and [010] 
axes) generate six types of AFM domains in MnTe films. 
Misfit strains in the MnTe/ZnTe multilayers reduce de- 



generacy to two possible domain types with the same q 
and mutually perpendicular L vectors in the film planei^ 

In a single-domain sample with the fixed magnetic or- 
der parameter (or AFM vector) L the equilibrium values 
of magnetostriction tensor, u, and electrical polarization, 
P, can be calculated by minimization of the following 
expression for the free energy 



F = / dV{[X'{uxx - Uyy) + 2\ifiUxy\{Li - Ly)0{z) + -Ui-jCijkiuu 



2eoK 



{PxUyz + PyU^x -I- PzU^y)} (7) 



r 



constructed from symmetry considerations. Parame- 
ters A are magnetoelastic constants related to the shear 
strain, c is a tensor of the elastic modula characteristic 
to cubic structure. Electrical properties are described by 
low-frequency dielectric constant Cpm = k + 1 and piezo- 
electric coefficient 614, ep is the vacuum permittivity. Pe- 
riodical alteration of the AFM (MnTe) and nonmagnetic 
(ZnTe) layers in a heterostructure is described by a form- 
function 9{z) which equals 1 in AFM layer and vanishes 
in a nonmagnetic spacer. 

In the case of ideal interfaces the strain induced by 
AFM ordering is homogeneous throughout the AFM and 
nonmagnetic layers and thus depends on the ratio of 
AFM layer thickness d to the period of superstructure 
D. In particular, macroscopic shear strain component 
which is responsible for polarization effect is given 



by the following expression 



u 



(spoil) 

xy ~ 



2c44 D ' 



(8) 



where p — ±1 distinguish between two orthogonal orien- 
tations of L||[100] and [010] in different domains, A^fp = 
)L] is sublattice magnetization. 

Strain-induced contribution'*^ into spontaneous polar- 
ization calculated from ([71), ([S]) is 



p(spoii) ^ _2e^^u 



(spon) 



euXiaW^p d 



C44 



D 



(9) 



For the multilayered structure the piezoelectric coeffi- 
cient eu{z) should be averaged over the AFM/NM bi- 
layer. 

It can be easily seen from the Expr. ^ that the di- 
rection of vector P correlates with that of AFM vector L 
(factor p) and takes an opposite direction in the domains 
with L||[100] and [010], in accordance with intuitive pre- 
dictions given above (see Fig. (S]). So, coupling between 
the ferroelectric and antiferromagnetic order is accom- 
plished by the rigid alignment of the antiferromagnetic 
axis perpendicular to the polarization direction. 

Characteristic value of the spontaneous ferroelectric 
polarization in a single-domain state can be evaluated 
from ([5]) and available experimental data (see Table WT\^ 
as p(>5p™) = 60 nC/cm2. For the thick MnTe film p(^P°") 
can be even higher because Mn atoms cause an en- 
hancement of electromechanical coupling in zinc-blende 
structures. "^^ It is remarkable that among the materials 
with comparable ferroelectric polarization (see Table |T| 
MnTe shows the highest Neel temperature. 

The value of the internal electric field £;(«p°'i) - 80 
kV/cm corresponding to the spontaneous polarization of 
MnTe is, in turn, close to strain-generated electric field 
in non-magnetic III-V heterostructures42i It seems to be 
large enough to be detected due to the change in the pho- 
toluminescence spectra, like it was done in the Ref. l20|. 
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V. SWITCHING OF POLARIZATION BY 
EXTERNAL MAGNETIC/ELECTRIC FIELD 

Practical applications of multifcrroics in information 
technology are based on their ability to maintain sin- 
gle ferroelectric domain state for a long time and change 
it under application of the external field. This can be 
achieved, for example, by using materials with a high 
ferroelectric Curie temperature and a robust large po- 
larization (say, RMnOs, R=Sc, Y, In, Ho-Lu, see Refs. 
141,42) that, in turn, gives rise to a large value of switch- 
ing field. In the systems with strong coupling between 
polarization and AFM order a single domain state can 
be easily fixed by a proper choice of the sample shape, 
even for small polarization. Moreover, if AFM and fer- 
roelectric domains are intimately related and match up 
spatially (as, e.g., should be in MnTe), electric polariza- 
tion of the sample can be switched by either magnetic or 
electric field. In this section we analyze field dependence 
of macroscopic polarization for different sample shapes. 

Existence of equilibrium domain structure in the fer- 
roics is usually attributed to the presence of long-range 
dipole-dipole interactions between the physically small 
ordered regions of the sample. Quantitative description 
of the phenomena is based on the account of shape- 
dependent contribution (stray energy) into free energy 
of the sample. In the case of ferroelectric (or ferromag- 
netic) materials this contribution can be written as 

^i-e='^{P,)N,k{Pk), (10) 
I 



where brackets (. . .) mean averaging over the sample vol- 
ume V and N is the second-rank depolarization tensor 
the components of which depend upon the shape of the 
sample. 

Equilibrium domain structure in AFM with nonzero 
magnetoelastic coupling can be described in a similar 
way. Long-range (destressing) effects originate from the 

local internal stress fields crj™'*^^ induced by magnetic 
ordering42iii Corresponding contribution into free en- 
ergy of the sample takes a form 

= I {(aj7^))H,,,,(4r'))}, (11) 

where i^jkpt is the forth-rank destressing tensor which, 
like N, depends upon the sample shape. 

It was already mentioned that the domain structure 
of MnTe/ZnTe multilayers is formed by two (out of six) 
types of the domains with (j) L||[100], P||[001], volume 
fraction a, and {ii) L|| [010], P|| [OOT], volume fraction (1— 
a). Switching between this two types can be performed 
by application of the external electric field E parallel to 
[001] axis or by the magnetic field H directed along one 
of the "easy" AFM axes ([100] or [010]). 

In the most practical applications^ macroscopic prop- 
erties (such as polarization or elongation) of the multido- 
main sample depend on the single parameter < a < 
1 that can be calculated from minimization of shape- 
dependent (<i>t_o + ^dostr) and field-dependent contribu- 
tions into free energy: 



$ = y 



Si - E,P., 



(spoil) 



a-- 



"-2 



(12) 



Here spontaneous polarization pj^P°"' within a domain 
is defined by the expression ([9]), x is the magnetic 
susceptibility of MnTe and 5*1 2 are shape-dependent 
coefficients^: 

51 = A'A/oV('"f)(Hii-H22), 

52 = 4Ar3 (pj'^P™))^ (13) 
+ Afo4[(A')'(Hii + H22 - 2H12) + 16A?eH66], 

where cr'™^' = ci™*'' = cr^™^^ are isotropic stresses in- 
duced in the film plane due to the lattice mismatch. 

Analysis of Eg. p^ shows that a single domain state 
of a sample (a = 1 or 1) can be achieved either by the 
proper choice of the shape (25*1 > S2) or application 
of the external field {E > Emu otc H > Hmb) where 
we've introduced characteristic monodomainization fields 



as follows: 

Emd = 7 \^ Hmd = \ ;-• (14) 

In the multidomain state macroscopic ferroelectric po- 
larization depends upon the external fields in a following 

way 

Y S2 Emd J^md J 

Let us analyze the properties of the coefficients 6*1.2 
assuming that i) the sample has the shape of an ellipse 
(with semiaxes a, b and eccentricity k = ^\ — b'^ /a?) 
within the film plane, film thickness is h; ii) elastic prop- 
erties of the material are isotropic (cn — C12 — 2C44) and 
are characterized with the shear modulus C44 and Poisson 
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ratio V = Cnj (cn + C12). The destressing and depolariz- 
ing coefficients can be calculated explicitly in two limiting 
cases: 

i) "pillar" with h ^ a > b: 
Si — 



S2 



2cmD (1 + VT^rfc2)2 

4A?6 - (A')2 



I -I' (1 + VT^P)4 

ii) "thin film" with h <ti b < a: 



(16) 



(17) 



where the dimensionless shape-factor is given by an 
integral 



J2(fc) = 



r''/^ (sin2 0-|-cos20/fc2)d0 

] \/l - sin^ (j) 

37rfcVl6, fc^O 
1, k-^l. 



(18) 



It is remarkable that coefficient ^2 that favors forma- 
tion of the domain structure is nonzero in both limit- 
ing cases, as can be seen from (fT3| . (fTB]) . (fT7|) . This 
fact can be easily extended to any geometry. Really, all 
the depolarizing effects (stray fields) are related with the 
flux of the corresponding (ferroelectric, magnetoelastic, 
etc) dipole moment through the sample surface. In our 
case the flux of ferroelectric polarization is nonzero only 
through the surface parallel to the film plane, thus "thin 
film" shows strong depolarizing ferroelectric effect while 
"pillar" shows none. In turn, AFM ordering produces 
stress dipoles within the film plane, so, corresponding 
flux is maximal through the side faces. As a result, strong 
destressing effect should be observed for the "pillar" , not 
for the "thin film" . For the intermediate case of the sam- 
ple with the developed side and face surfaces ("ball" or 
"cube" ) coefficient S2 is contributed by both depolarizing 
mechanisms and so is nonvanishing for any shape. 

Calculations based on the formulas p^ . p?)) and ex- 
perimental data from Table IVIII show that in the case 
of MnTe multiferroic the destressing effects are much 
stronger (^2 ~ 10^ J/m^) that the electric depolariza- 
tion (^2 ~ 5 • 10^ J/m"^), so, monodomainization can be 
easier achieved in the "thin film" . 

Another interesting feature of the destressing phenom- 
ena is existence of the effective internal shape-induced 
field described by coefficient Si. It arises due to cross- 
correlation between isotropic (induced either by sub- 
strate or by magnetovolume effect) and anisotropic in- 
ternal stresses and has no analog in ferromagnetic or fer- 
roelectric materials. Coefficient Si, as seen from (IT| 



dnl), (Ull), depends upon the eccentricity k and vanishes 
for the samples isotropic within the film plane (a — b). 

So, we can deduce that the shape of the sample influ- 
ences the domain structure in two ways. First, the field of 
monodomainization (jl4p strongly depends upon the h/b 
ratio. Characteristic values of the electric and magnetic 
monodomainization fields for two limiting cases {h ^ b 
and h <^ b) are given in Table IVnl The value of Hmd 
for the thick sample ("pillar") is close to the value of 
anisotropy field 3.25 T (see Reflisl) and diminishes down 
to 0.8 T for 'thin film". 

Second, single domain state is energetically favorable 
for the samples with the overcritical eccentricity (calcu- 
lated from the condition 25*1 > 5'2)- For example, for the 
thick sample the critical eccentricity is 0.57 (correspond- 
ing aspect ratio a/b = 1.22). 

Shape effects are clearly seen in Fig. IH] which shows 
field dependence of macroscopic ferroelectric polarization 
for the "thin film" with h/b = 0.1 and for the "piUar", 
both samples having the same in-plane eccentricity 0.1 
(corresponding aspect ratio a/b = 1.005) In the first case 
monodomainization field is smaller and biasing effect is 
more pronounced. 

Monodomainization of the sample can be also achieved 
by the combined application of the electric and magnetic 
field, as illustrated in Fig. [7] for the case of "thin film" 
with zero eccentricity. Relation between the switching 
fields in this case is given by the formula 



Hy - 



Emd 



tt2 



= 1. 



(19) 
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N e 




6 4 2 


-20- 
-40- 






-80- 





Figure 6: (Color online) Macroscopic ferroelectric polariza- 
tion as a function of external magnetic field calculated ac- 
cording to (|15l) for "thin film" with h/b = 0.1(curve A, ma- 
genta) and "pillar" (curve B, dark cyan). The eccentricity in 
both cases is fc = 0.1. Magnetic field H is switched between x 
and y directions. Dotted curve (magenta) shows the possible 
field switching between two opposite polarizations in a single 
domain sample. 
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Figure 7: (Color online) Phase diagram in E — H variables for 
the "thin film" {h/b — 0.1) with zero eccentricity. Magnetic 
field H is switched between x and y directions. 



VI. VISUALIZATION OF THE DOMAIN 
STRUCTURE 



The zinc-blend wide-gap semiconductors are also 
known to show nonlinear optical properties studied in 
second harmonic generation (SHG) experiments."*^ In the 
paramagnetic phase these crystals have the only nonzero 
SHG tensor component (ii4(— 2cj, cj, w) — ^25 = dse 
which is rather large, e.g. for ZnTc it is 119 pm/V at 
the fundamental wavelength 1047 nm."*^ AFM ordering 
may bring to being additional SHG components. To cal- 
culate them we analyze the response to the external elec- 
tric field with account of the nonlinear contribution into 
free energy 



p(nonlin) _ _^^p p p , j/_/p2 _ p2w _ ) + f P P „ 4- P P u +PPu ) 

^ — 23^ y '^ \ X y ^xx '-'"yy / ' x-^ y ^xy \ y-^ z ^yz \ ^ z^ x^zx J 

+ e{z)[^,{Ml + M2)P,(P2 „ p2) ^ ^^(^,^2 _ ^^2)p3]^ (20) 



Coefficients 77', 7744 introduced in ([^ describe elec- 
trostrictive effect and the last term describes high-order 
magneto-electric coupling. 

Calculations show that in the AFM phase the crystal 
becomes biaxial with anisotropic dielectric tensor: 

e,, = ^M2p(ep„.-1), (21) 

C44 

in accordance to symmetry prediction for Ci point group 
of AFM phase. The SHG tensor components propor- 
tional to AFM vector are of two types. Like dielectric 
coefficients, the components 

4, + = 47744rfi4Ai6 2^ ^ 24^P°")r;44di4P, (22) 
C44 " 

have opposite signs for the domains with perpendicular 
orientation of AFM vectors (p = il)- The other coeffi- 
cients are insensitive to the domain structure: 

dzz = -Mo'6eo(epm(c^) - 1). (23) 

The symmetry predicted difference between the (ii4, (^25 
and (iae components is due to the high-order terms in the 
magnetic order parameter and is neglected. 

Rotation of the polarization plane which may stem 
from anisotropy of dielectric tensor and effect of SHG are 
sensitive to the direction of AFM vector and thus open 
a possibility to visualize the AFM domain patterni^ 



VII. SUMMARY AND CONCLUSIONS 

In the present paper we demonstrate for the first time 
that MnTe epitaxial films and MnTe/ZnTe multilayers 
grown in a convenient (OOl)-direction should become 
electrically polarized below the Neel temperature. Spon- 
taneous ferroelectric polarization results from the bend- 
ing of highly ionic Mn-Te-Mn bonds induced by magne- 
tostrictive shear strain. 

In contrast to many other multiferroics, MnTe has no 
macroscopic magnetic moment and the direction of fer- 
roelectric polarization is coupled with the orientation of 
AFM order parameter. Thus, the domain structure is 
formed by the combining action of the electric and mag- 
netoelastic dipole-dipole interactions. Corresponding de- 
polarizing long-range contribution into free energy of the 
sample is nonzero for any sample shape. In the case 
of thick film depolarizing effects are governed mainly 
by magnetoelastic mechanism and corresponding mon- 
odomainization field is of the same order as characteristic 
anisotropy field. In the case of thin film depolarization is 
due to the presence of the electric dipole interactions and 
corresponding monodomainization field is much smaller. 

Due to the clamping between the electrical and AFM 
domains switching between the opposite direction of fer- 
roelectric polarization can be induced by application of 
the magnetic field parallel to "easy" AFM axis. Vice 
versa, one can switch between different (lOO-oriented and 
[010]-oriented AFM domains) by application of the elec- 
tric field. This effect seems to be useful for controlling 
the state of the adjacent ferromagnetic layer coupled to 
the multiferroic through exchange interactions like it was 
proposed in Refi9|. 
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AFM ordering may also induce additional (anisotropy) 
components of dielectric and SHG tensors. The sign of 
the effect is sensitive to orientation of AFM vectors and 
thus opens a possibility to distinguish different domain 
types. The spontaneous ferroelectric effect is peculiar to 
the collinear AFM- III ordering and is not allowed (from 
symmetry point of view) in the canted AFM-III struc- 
ture (e.g., such an effect is impossible in MnTe/CdTe 
multilayers). 
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Table III: Values of constants for MnTe (ZB) used in calcula- 
tions and calculated (see text for notations). 



Constant 


Value 


Source /formula 


Rom. 


Tn 


65 K 


Ref.23 




c/a 


1.06 


Ref.23 






0,003 


Ref.23 




ei4 


0.097 C/m^ 


Ref.39 


CdgiMnigTe 


^pm 


9,5 


Ref.39 




Cll 


53.6 GPa 


Ref.39 


CdgiMmgTe 


Cl2 


37 GPa 


Ref.39 


CdsiMnigTe 


C44 


18.2 GPa 


Ref.39 


CdsiMmgTe 


Mo 


615 Gs 


Ref.21 




X 


1.8-10-* cgs 


Ref.47 


r = 


V 


0.41 


Cl2/(cil + C12) 




Un-ii 


0.03 


0.5(c/a - 1) 




0"mf 


2.7 GPa 


(Cll + Cl2)ttmf 




X'Mg 


0.06 GPa 


C44(^Uxx '^yy) 




p(spon) 
^ z 


60 nC/cm^ 


® 




^(spon) 


80 kV/cm 


Pi=^°"V(eo«) 




-Emd 


750 kV/cm 


m 


"pillar" 


-Emd 


40 kV/cm 


m 


"thin film" 




3.5 T 


m 


"pillar" 




0.8 T 


m 


"thin film" 
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anisotropic exchange mechanism in this system. 
'^^ The analogous contribution into polarization may also 
arise from the direct magneto-electric coupling omitted in 
the Expr.Q. 

For the ellipsoid-shaped samples disregarding small con- 
tribution from the domain walls and possible closure do- 
mains.) 

We use Voight notations for the components of the sym- 
metric 2-nd and 4-th rank tensors 



